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INTRODUCTION 

This  report  presents  a  method  for  estimating  the  coordinates  of  a 
moving  target  as  a  function  of  bearing  direction  cosines  measured  from 
a  tracking  platform.  The  procedure  can  be  considered  as  a  generaliza¬ 
tion  of  Ekelund  ranging  since  parameter  estimates  are  computed  after 
bearing  measurements  are  taken  during  a  period  when  the  tracker  maneuvers. 
The  key  assumption  for  the  method  is  that  the  target  is  moving  at  a  con¬ 
stant  speed  on  a  fixed  heading.  This  is  the  same  assumption  made  for 
Ekelund  ranging.  In  contrast,  however,  the  method  which  is  the  subject 
of  this  report  does  not  constrain  the  course  of  the  tracking  platform. 


DEVELOPMENT  OF  THE  TRACKING  PARAMETERS 


For  a  fixed  coordinate  system,  let  x^,(t)  and  x^(t)  denote  the  x  coordi¬ 
nates  of  the  target  and  tracker  at  time  t.  Let  y^,(t)  and  Yg(t)  denote  their 
y  coordinates.  Let  B(t)  denote  the  true  target  bearing  at  time  t  measured 
with  respect  to  the  y  axis.  For  example,  a  target  on  the  y  axis  would  be 
at  0  or  180  degrees.  Then,  from  Figure  1, 


and 


sin  B(t)  = 


xT(t)  -  xfi(t) 
R(t) 


cos  B(t)  = 


yT(t)  "  yB(t) 
R(tj 


(la) 


(lb) 


where  R(t)  is  the  true  target  range. 
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Assume  that  the  tracker  uses  a  circular  hydrophone  array  to  detect 
acoustic  waves  radiating  from  the  target.  Bearing  information  is  usually 
obtained  by  delay-and-sum  beamforming.  Remember  that  the  direction  cosines 
sin  B  and  cos  B  are  needed  to  determine  the  delays  used  to  steer  a  beam  in 
direction  B.  Suppose  that  the  beam  angle  that  gives  maximum  signal  energy 
during  the  integration  time  is  B.  This  is  the  estimate  of  the  target  bear¬ 
ing  during  the  integration  time.  Let  s(t)  -  sin  B(t),  c(t)  =  cos  B(t), 
and  s(t)  and  c(t)  denote  the  direction  cosines  corresponding  to  B.  They 
can  be  obtained  from  the  signal  processor  of  the  beamformer.  Levin1  and 
and  Hinich  and  Shaman2  show  that  these  estimators  of  sin  B  and  cos  B  are 
maximum  likelihood  estimators  if  the  ambient  noise  is  Gaussian  and  spatially 
incoherent.  If  the  array  gain  is  large,  moreover,  these  estimators  are 
approximately  Gaussian,  unbiased,  and  independent.  For  a  circular  array 
geometry,  the  variances  of  s(t)  and  c(t)  are  equal  and  are  inversely  pro¬ 
portional  to  the  energy  signal-to-noise  ratio  (SNR).  Expressions  for 
the  bearing  and  direction  cosine  errors  as  a  function  of  SNR,  aperture, 
and  the  number  of  hydrophones  is  given  by  McDonald  and  Schultheiss ,3 
Clay,  Hinich,  and  Shaman,4  and  Hinich.3 

Suppose  that  the  tracker  estimates  sin  B(t)  and  cos  B(t)  at  discrete 
time  points  t^  =  nx  where  I  is  the  integration  time  of  the  beamformer. 


*M.  J.  Levin,  "Least-Squares  Array  Processing  for  Signals  of  Unknown  Forms," 
Radio  and  Electronic  Engineer,  Vol.  29,  pp.  213-222  (1965). 

2M.  J.  Hinich  and  P.  Shaman,  "Parameter  Estimation  for  an  R-Dimensional 
Plane  Wave  Observed  with  Additive  Independent  Gaussian  Errors,"  The  Annals 
of  Mathematical  Statistics,  Vol.  43,  pp.  153-169  (1972). 

3V.  H.  MacDonald  and  P.  M.  Schulthe  iss,  "Optimum  Passive  Bearing  Estimation," 
Journal  of  the  Acoustical  Society  of  America,  Vol.  46,  pp.  37-43  (1969). 

4C.  S.  Clay,  M.  J.  Hinich,  and  P.  Shaman,  "Error  Analysis  of  Velocity  and 
Direction  Measurements  of  Plane  Waves  Using  Thick  Large-Aperture  Arrays," 
Journal  of  the  Acoustical  Society  of  America,  Vol.  53,  pp.  1161-1166, 

(1973). 

5M.  J.  Hinich,  "On  Errors  in  Some  Papers  on  Array  Processing,"  Journal  of 
the  Acoustical  Society  of  America,  Vol  65,  pp .  530-531  (1979). 
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The  statistical  expressions  derived  in  the  next  section  are  simplified 
if  we  set  the  time  origin  in  the  middle  of  the  sampling  period;  i.e.,  let 
n  assume  integer  values  n  =  -N/2,  -N/2  +  l,...,N/2  -  1,  N/2  where  N  is 
even.  Thus,  the  sampling  interval  is  (N  +  1)  T.  Assume  that  T  is  selected 
to  ensure  that  the  estimates  are  uncorrelated  over  time. 

Several  important  assumptions  will  now  be  made  about  the  target  motion 
and  the  SNR  during  the  sampling  period.  First,  assume  that  the  target's 
velocity  v^  and  its  heading  are  constant.  Thus,  for  -N/2  <  n  <  N/2 

x_  (t  )  =  xT  (0)  +  v  t  sin  a_  (2a) 

Ini  In  I 

and 

yT  =  yT  ^  +  VT  tn  C0S  ffT  '  (2b) 

Second,  assume  that  vjt  <<  R(0)  and  v  NT  <<  R(0),  where  v  is  the  average 

Id  d 

speed  of  the  tracker.  This  implies  that  the  range  R(t)  is  approximately 
constant  during  the  sampling  period.  Now  let  R  denote  the  average  range 
in  the  sampling  period.  Finally,  assume  that  the  SNR  varies  sufficiently 
slowly  during  this  period  so  that  the  SNR  can  reasonably  be  approximated 
by  a  constant.  This  assumption  implies  that  the  variances  of  the  direc¬ 
tion  cosines  are  approximately  constant.  This  assumption  will  be  relaxed 
in  a  later  section. 

To  simplify  notation,  select  the  time  unit  so  that  T  =  1  and  thus 
tfi  =  n.  It  then  follows  from  Equations  (1)  and  (2)  and  the  above  assump¬ 
tions  that 

s(n)  =  R  1  xT(0)  +  (R  1  vTsin  (*T)n  -  R  1  xg(n)  +  tg(n)  (3a) 

and 
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c(n)  =  R  1  yT(0)  +  (R  1  vTcos  aT)n  -  R  1  yB(n)  +  ec(n)  (3b) 

where  the  errors  have  the  following  properties  (for  circular  arrays  whose 
gain  is  large): 

1.  £g(n)  and  6  (n)  are  independent  Gaussian  random  variables 
with  a  common  variance  denoted  a2 

2.  es(n)  an<*  £g(n')  (and  Ec(n)  and  £c(n'))  are  uncorrelated 
for  all  n  ^  n' 

3.  The  expected  values  of  these  errors  are  approximately  zero 
if  R(t)  £  R  during  the  sampling  period. 

Now,  s(n) ,  c(n),  xfi(n)  and  yfi(n)  are  observed  for  n  =  -N/2, . . . ,N/2. 
Hence,  the  coordinates  [x^,(no),  y^,(nQ)]  of  the  target  at  time  nQ  can  be 
estimated. 

LEAST  SQUARES  ESTIMATES  OF  TRANSFORMED  PARAMETERS 

The  maximum  likelihood  estimation  of  these  target  parameters  is 
facilitated  by  the  following  transformations: 


a  =  R_1 

X 

xT(0) 

a  =  R*1 

y 

yT(o) 

„  -l 

„  -1 

p  =  R 
rx 

vT  sin  aT 

p  =  R 

y 

vT  cos 

and 

b  =  -R'1  .  (5) 


5 


Thus,  from  Equation  (3) 


s(n)  =  n  +  b  x^(n)  +  £s(n)  (6a) 

and 

c(n)  =  ay  +  Pyn  +  b  yB(n)  +  ec(n)  (6b) 


for  n  =  -N/2, . . . ,N/2.  Since  the  errors  are  Gaussian  and  have  constant 
variance,  the  ordinary  least-squares  (OLS)  estimators  of  a^,  a^,  p^,  and 
py  are  maximum  likelihood.6  They  also  have  a  joint  Gaussian  distribution. 
The  maximum  likelihood  estimator  of  b  is  a  weighted  average  of  b^  and  b^, 
the  OLS  estimators  of  b  as  computed  from  Equations  (6a)  and  (6b).  The  maxi¬ 
mum  likelihood  estimators  of  the  target  coordinates  are  functions  of  these 
OLS  estimators  and  are  computed  from  the  transforms  defined  in  Equation  (4). 
These  estimators  are  presented  in  the  next  section,  Target  Range  and 
Coordinate  Estimates.  Before  presenting  the  weights  needed  to  compute 

A  A  A 

b,  the  OLS  estimators  a^,  p^,  and  b^  and  their  statistical  properties 
are  presented. 


The  expressions  for  the  OLS  estimators  of  a^,  P^,  and  b^  are  simplified 
if  the  origin  of  the  coordinate  system  is  placed  at  the  centroid  of  the  plat¬ 
form's  track  during  the  sampling  period.  If  this  is  done,  then  <Lx_(n)  = 

D 

ly  (n)  =  0.  The  OLS  estimator  of  a  is 
l5  X 

5x  =  (N  +  l)"1  Is(n)  .  (7) 

Its  variance  is  simply 

o2  =  (N  +  l)"1  o2  .  (8) 

a 

x 


6M.  G.  Kendall,  The  Advanced  Theory  of  Statistics,  Chapter  22,  New  York: 
Hafner,  Third  Edition  (1951). 
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Now  define  s(n)  =  s(n)  -  s,  where  s  =  (N  +  1)  Is(n),  the  mean  of 
s(n).  The  OLS  estimators  8  and  b  can  be  expressed  in  vector  and  matrix 

X  X 

fo>-m  as  follows: 


In2  Inx_(.n) 

D 


b  J  'InxB  IxB(n) 


Ins (n) 


IxB(n)  s(n) 


The  sums  are  taken  from  n  =  -  N/2,...,N/2.  Thus, 


8X  =  Dx  [lx2 (n)  Ins (n)  -  Inxg(n)  Ixg(n)  s(n)] 


b  =  D  1  [In2  lx  (n)  s(n)  -  InxD(n)  Ins(n)] 

X  X  D  D 


where 


=  In2  Ix2(n)  -  (InXg(n) ] 2 


From  the  triangle  inequality,  D  =0  if,  and  only  if,  x_(n)  is  a  linear 

X  D 

function  of  n  -  t  .  Thus,  d^  f  0  if  the  tracker  changes  course  or  speed 
during  the  sampling  period. 


The  variance-covariance  matrix  of  B  and  b  is 

rx  x 


Ix2(n) 


-InXg(n) 


■InXg(n) 
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Thus,  the  variances  of  8  and  are 

x 


°8  =  °2  Dx 

px 


and 


a 2  =  a2  D  1  In2 

b  x 

x 


It  can  be  shown  that  a  is  uncorrelated  with  B  and  b  . 

X  XX 


(U) 


(15) 


It  should  be  obvious  from  Equations  (6a)  and  (6b)  that  the  OLS  estimators 

A  A  A  A 

a  ,  B  i  and  b  are  similar  to  a  ,  B  >  and  b  with  y  (n)  and  c(n)  in  place 

place  of  x  (n)  and  s(n).  Thus  a  =  (N  +  1)  ^  Ic( n) . 

B  y 

By  =  Dy1  [^Yg(n)  Inc(n)  -  Inyg(n)  IyB(n)  c(n)],  (16) 

and 

6^  =  Dy1  [In2  IyB(n)c(n)  -  Inyfi(n)Inc (n) ] ,  (17) 


where  c(n)  =  c(n)  -  c  and 

Dy  =  In2  ly2  (n)  -  (InyB(n)]2  .  (18) 

Again,  Dy  t  0  if  the  tracker  makes  a  course  or  speed  change. 

The  estimators  a  ,  B  >  and  b  are  independent  of  a  ,  B  >  and  b 

xxx  y  y  y 

since  the  errors  (e^Cn)]  are  independent  of  (e  (n)].  It  then  follows 
from  statistical  theory6  that  b,  the  maximum  likelib  od  estimator  of  b, 
is  the  weighted  average 


*  O'1 


b 

x 


6  ibid. 
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TARGET  RANGE  AND  COORDINATE  ESTIMATES 

Except  in  some  special  cases,  the  variances  of  the  OLS  estimators 
go  to  zero  as  N  ■>  ®.  The  rate  at  which  this  occurs  is  a  function  of  N 
and  depends  upon  the  form  of  [xg(n),  yfi(n) ] ;  i.e.,  the  form  of  the  plat¬ 
form's  track.  The  variances  also  go  to  zero  as  o'c ,  the  variance  of  the 
bearing  direction  cosine  errors,  goes  to  zero.  As  is  shown  in  the  follow- 
theorem,  the  maximum  likelihood  estimate  of  the  range  is  biased  if  Ra^  is 
not  small.  Assume  then  that  N  is  sufficiently  large  and  a  is  sufficiently 

small  so  that  Ra.  is  small, 
b 

Theorem  1.  R  =  -1/b  is  the  maximum  likelihood  estimator  of  R.  The 

A  A 

bias  in  R  due  to  the  nonlinear  transformation  of  b  is 

E(R-R)  =  R3cj2  +  0  (R3a£)  (21) 
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and  its  approximate  root  mean  square  error  is 

rmse  (R)  =  R2o.  +  0  (R3o?),  (22) 

D  b 

where  o2  is  given  by  Equation  (20).  Thus,  the  bias  is  an  order  of  magni¬ 
tude  smaller  than  the  rmse  when  Ro,  <<  1. 

D 

Proof,  b  =  b  +  £.  is  the  maximum  likelihood  estimator  of  b  =  -R 
-  b  ~ 

If  f(b)  is  a  continuously  different  ble  function  of  b,  then  f(b)  is  the 
maximum  likelihood  estimator  of  f(b).  Consequently,  R  is  the  maximum 
likelihood  estimator  of  R. 

The  error  is  Gaussian  N(0,o2).  Thus  Ee^  =  Ee3  =  0,  Ee2 
Ee^  =  3oJ.  Since 

R  _  1 

R  1  -  Reb 

=  1  +  Reb  +  R2e2  +  Raeb  +  . . .  , 

it  follows  that  the  expected  value  of  R/R  is 
E(R/R)  =  1  +  R2afe  +  0  l(Rob)4], 
and  its  mean  square  error  is 

mse  (R/R)  =  R2o2  +  0  [(Rab)4]. 

Expressions  (21)  and  (22)  follow  from  Equations  (24)  and  (25). 

Given  R,  the  maximum  likelihood  estimators  of  x_(n  )  and  y,_(n  ) 

l  o  l  o 

can  easily  be  obtained.  Accurate  estimates  of  these  target  coordinates 
require  that  Rab  be  small. 

Theorem  2.  The  maximum  likelihood  estimators  of  [x^,(no),  y^. (nQ ) ] 

are 


=  Ob ,  and 

(23) 

(24) 

(25) 
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iT(no)  =  R(Ax  +  “o^ 


and 


(26) 


yT(n  )  =  R(a  +  n  p  ) 
T  o  y  ory 


They  are  independent.  For  large  N,  their  distributions  are  approximately 
Gaussian  with  zero  means  and  root  mean  square  errors 


rmse  ($T)  £  Ro[J  +  x2^) 


In2_ 
D 

x  y 


Inx_(n)  2x2(n)  . 

-  W.)  r-hr  * it—  i1/2 

x  y  x 


and 


rmse  (yT)  -  Ro(i  +  y2(no) 

x  y 

In yEU)  2  2y|(n)  1/2 

2noyT<‘no^  D  +  D  nr>  D 

x  y  y 

Proof:  Write  a  =  a  +  e  ,P  =  p  +e„,  and  b  =  b  +  e.  . 

x  x  a  ’  *x  Kx  p  ’  b 

x  x 

e  and  e„  are  uneorrelated  Gaussian  errors.  The  error  e,  is  correlated 
a  p  b 

x  x 

with  e„  .  Apply  Equations  (4)  and  (23)  to  Equation  (26)  to  obtain  the 
px 

following  approximation  (in  the  errors): 

xT(nQ)  ~  R(1  +  Reb)R  1xTCO)  +  ea  +  nQ(R  ^sin  aT  +  e„  ) 

x  px 

=  (1  +  Reb)(xT(nQ)  +  Refl  +  Rn^g  ] 

x  px 


£  xT(n  )  +  R(e  +  xT(n  )e,  +  n  e„  ) 
^  T  o  a  Tod  op 

x  rx 


(27) 
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Apply  Equations  (14),  (20),  and  the  expression  for  the  covariance  between 
b  and  to  Equation  (27)  to  derive  the  expression  for  rmse  (x^,).  Derive 
the  expression  for  rmse  (y^)  in  a  similar  manner. 


NON-CONSTANT  BEARING  VARIANCE 


Now  relax  the  restriction  that  the  SNR  is  constant  during  the  sampling 
period.  There  are  several  methods  for  estimating  the  SNR  associated  with 
a  specific  bearing  estimate  B.  One  method  uses  output  from  the  beamformer. 
If  the  noise  field  is  isotropic  in  the  sector  B-6<B<B+6,  then  the 
average  of  the  energy  in  the  beams  spanning  this  26  wide  sector  (excluding 

a* 

B)  is  an  estimate  of  the  noise  field  energy.  Let  e(6)  denote  this  esti- 

A 

mate.  The  energy  in  the  B  beam  is  an  estimate  of  signal  plus  noise.  Let 
e(B)  denote  this  estimate.  Then 


P 


_  e(B) 
"  e(6) 


1 


(28) 


A  somewhat  more  precise  estimate  uses  the  average  of  the  coherence 
between  pairs  of  sensor  channels.  To  illustrate  this  method,  let  y2 
denote  the  average  estimated  square  coherence  between  two  hydrophone 
channels  over  the  bandwidth  of  the  signal.  Then  the  SNR  is  estimated 
by 


The  maximum  likelihood  estimators  are  functions  of  the  variance  of 
the  direction  cosiqes  if  these  variances  change  during  the  sampling 
period.  Since  the  variance  depends  on  the  SNR,  it  must  be  estimated  every 
time  a  bearing  is  obtained  if  one  desires  to  approximate  the  statistical 
properties  of  maximum  likelihood  estimators. 
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The  adjustment  for  a  changing  error  variance  is  easy  using  the  linear 
model  approach;  i.e.,  use  weighted  least  squares.7  Simply  divide  the  inde¬ 
pendent  variables  [n,  x_(n),  and  y  (n)J  and  the  dependent  variables  |s(n) 

D  D 

and  c(n)|  by  the  estimated  variance  for  each  n.  Since  the  variance  is 
inversely  proportional  to  the  SNR  given  a  narrowband  signal,  the  adjustment 
is  made  by  multiplying  the  variables  by  p(n) ,  a  consistent  estimator  of  the 
SNR  at  time  t  =  nt .  All  the  expressions  for  the  estimators  and  their 
properties  given  in  the  previous  sections  hold  if  n,  xD(n),  y_(n),  s(n), 

O  D 

and  c(n)  are  multiplied  by  p(n).  If  the  true  SNR  for  each  n  were  known, 
then  the  adjusted  estimators  would  be  maximum  likelihood. 

It  is  important  to  remember  that  p  is  an  estimate  of  the  SNR  and  not 
the  true  SNR  for  any  type  of  estimator.  The  error  inherent  in  p  increases 
the  error  m  the  parameter  estimates.  Consequently,  it  is  generally  better 
to  use  the  oilini'”  least  squares  approach  rather  than  weighted  least  squares 
(the  multiplication  adjutment)  when  the  SNR  is  slowly  varying  during  the 
sampling  period. 


The  following  ad  hoc  compromise  between  ordinary  least  squares  and 
weighted  least  squares  estimation  may  provide  a  more  robust  estimation 
method.  The  compromise  rejects  a  bearing  estimate  if  its  estimated  SNR 
is  below  some  threshold  value.  Suppose,  for  example,  that  the  SNR  for 
s(2)  and  c(2)  is  below  the  threshold.  Then  delete  s(2),  c(2),  x  (2), 

»  n 

y_(2),  and  2  from  the  sums  in  the  expressions  for  6  ,  B  ,  b  ,  and  b  . 

r  rx  y  x  y 

The  target  parameters  are  estimated  from  those  bearings  whose  SNR  exceed 
the  threshold,  using  the  ordinary  least  squares  approach. 


SIMULATION 

It  is  assumed  the  target  is  moving  at  a  constant  velocity  of  90  yards/ 
15  seconds  or  0.36  kyd/minute  on  a  fixed  heading  of  90  degrees.  For 


7f.  A.  Graybili,  An  Introduction  to  Linear  Statistical  Models,  Vol .  1, 
New  York:  John  Wiley  &  Sons,  Inc.  (1961). 
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convenience,  the  tracker's  speed  is  set  at  the  same  rate.  This  is  not 
necessary  for  the  results  to  hold;  it  merely  facilitates  the  calculations. 

Upon  observing  the  target,  the  tracker  begins  its  tracking  maneuver. 

In  the  simulation,  the  track  is  circular  (Figure  2).  The  radius  of  the 
circle  is  such  that  the  tracker  completes  two  circles  during  the  tracking 
sequence  at  the  constant  velocity  of  0.36  kyd/minute.  This  particular 
track  was  chosen  for  several  reasons.  First,  it  is  symmetric;  that  is, 
one  circle  is  completed  in  each  half  of  the  cracking  sequence.  Thus, 
the  centroid  of  the  tracking  sequence  is  the  origin  of  the  coordinate 
system.  Second,  the  track  of  two  circles  yields  two  estimates  of  the 
target's  bearing  and,  hence,  the  estimated  bearing  sine  and  osine,  at 
each  point  where  the  tracker  takes  a  reading.  This  enables  the  OLS 
routine  to  estimate  more  accurately  the  target's  heading  since  any  differ¬ 
ence  in  the  two  estimates  of  the  target's  bearing  is  due  to  a  change  in 
the  target's  position.  And  finally,  a  circular  track  minimizes  the  expected 
variance  of  the  target's  range  for  a  smooth  track  when  the  tracker  has  no 
prior  knowledge  of  the  target's  heading  In  this  simulation,  45  observa¬ 
tions  of  the  target's  estimated  bearing  are  taken.  This  number  is  strictly 
arbitrary  although  the  routine  requires  that  the  number  of  observations  be 
an  odd  integer.  Hence,  the  radius  of  the  circular  track  is  45 /Un  or  0.322 
yard.  The  tracker,  therefore,  travels  3.96  kyds  after  having  first  observed 
the  target  taking  an  observation  every  15  seconds  or  every  90  yards. 


Each  estimate  of  the  target's  bearing  is  used  to  calculate  an  estimate 
of  s(n)  and  c(n).  The  resulting  45  estimates  are  used  to  estimate  the 
parameters  of  Equations  (6a)  and  (6b)  by  means  of  OLS.  The  estimated 
parameters  obtained  are  used  to  calculate  estimates  of  the  target's  head¬ 
ing,  range,  and  velocity.  The  estimated  heading  is  d_  =  n/2  -  tan  *  (p  /p  ) 

i  y  x 

It  is  computed  in  this  manner  because  P^  is  close  to  zero,  and  thus  small 
variations  in  P  produce  large  fluctuations  in  the  estimate  tan  ^(P  /P  ). 

y  *  y 

The  estimated  range,  R,  is  calculated  according  to  Equation  (5).  However, 

early  work  with  the  simulation  suggested  that  use  of  weighted  estimator  b 

provided  highly  variable  and  unrealistic  estimates  of  the  range.  This  is 

due  to  the  fact  that  b  ,  like  P  ,  is  rlose  to  zero  and  is  often  estimated 

y  y 
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as  a  positive  fraction.  The  smaller  of  b  and  b  is  selected  to  estimate 

x  y 

the  range.  The  program  consistently  selected  b^  since  is  90  degrees. 
The  weighted  estimator  b  is  only  optimal  when  the  sample  size  is  large. 


The  sample  size  is  not  large  in  this  simulation.  The  results  in  Table  1 

~-l 

are  derived  by  using  the  weighted  estimator  of  range,  -b  .  Table  2 

~-l 

presents  the  results  obtained  by  using  the  ad  hoc  estimator,  -b  The 

A  A  >.  I/O  ^ 

estimated  velocity  v^,  is  estimated  by  v^  =  [R2(p2  +  p2)] 


These  estimates  are  then  used  to  compute  the  estimated  target  coordi¬ 
nates  at  the  end  of  the  tracking  sequence  x^,(N/2)  and  y^,(N/2)  according  to 
the  formula  of  Theorem  2.  In  addition,  the  theoretical  rmse's  of  these 
values  are  calculated.  They  are  denoted  rmse^,(x^)  and  rmse^,(y^,)  in  the 
tables  and  are  presented  for  comparison  to  the  actual  values  found  in  the 
simulation. 


The  simulation  is  repeated  400  times.  The  tables  present  the  actual 
values  as  well  as  the  mean,  mean  bias,  standard  deviation,  and  the  root 
mean  square  error  (rmse)  for  each  of  the  estimated  variables.  The  estimated 
variables  are  (1)  the  target's  estimated  coordinates  at  the  end  of  the 
tracking  sequence,  x^,(N/2)  and  y^,(N/2,  (2)  the  target's  estimated  average 
range  during  tracking,  R,  (3)  the  target’s  estimated  heading,  or^. ;  and  (4)  the 
target's  estimated  velocity,  v^.  In  addition,  the  actual  maximum  and  mini¬ 
mum  range  to  the  target  during  tracking  is  presented. 


Compare  the  estimates  in  Table  2,  which  are  obtained  using  the  ad 
hoc  estimator  of  range,  -b  ^ ,  with  those  of  Table  1,  which  were  obtained 
using  the  weighted  estimator  of  range.  It  is  apparent  that  the  use  of 
the  weighted  estimator  significantly  increases  both  the  bias  and  variance 


of  all  the  estimates  obtained.  The  exception  is  the  estimated  heading 


a^,.  This  result  follows  from  the  fact  that  the  weighted  estimator  con¬ 
sistently  over  estimates  the  range  which  is  used  to  calculate  all  estimates 
except  that  of  the  estimated  heading  a^,. 
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TABLE  1 


RESULTS  OBTAINED  USING  WEIGHTED  ESTIMATOR  OF  RANGE  o  =0.20 


End  of  Track 

Actual 

Xj.(N/2) 

1.9799986 

yT(N/2) 

19.7999878 

I ^(N/2) 

19.8363953 

cL 

T 

90 . 0000000 

VT 

0 . 3600000 

Mean 

Mean  Bias 

3.6655350 

1.6855354 

36.7997894 

16.9997864 

36.8514709 

17.0150604 

89.7862091 

-0.2137576 

0.6682853 

0.3082853 

e 


Standard 

Deviation 

RMSE 

0.2367992 

1.7020874 

0.1459860 

17.0004120 

2.3269196 

17.1734467 

0.5106176 

0.5535945 

0.0427839 

0.3112399 

Maximum  and  minimum  range  during  tracking:  20.181183  and  19.485062 


N  =  45  R  =  19.80  rmse^,(xT)  =  0.069  rmseT(yT)  =  0.667 


TABLE  2 

RESULTS  OBTAINED  USING  AD  HOC  ESTIMATOR  OF  RANGE  o  =0.20 

£ 

Standard 


End  of  Track 

Actual 

Mean 

Mean  Bias 

Deviation 

RMSE 

*t(N/2) 

1.9799986 

1.9686899 

-0.0113090 

0.0916492 

0.0923442 

yT(N/2) 

19.7999878 

19.7851410 

-0.0148508 

0.0908208 

0.0920270 

Rj(N/2) 

19.8363953 

19.8159637 

-0.0204318 

0.8945770 

0.8948103 

®T 

90 . 0000000 

89.8274689 

-0.1725032 

0.5489491 

0.5754150 

0 . 3600000 

0.3588170 

-0.0011830 

0.0168327 

0.0168743 

Maximum  and  minimum  range  during  tracking:  20.181183  and  19.485062 
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Table  3  shows  the  results  when  the  bearing  sine  and  cosine  errors  are 
doubled  from  0.2  degree  in  Table  2  to  0.4  degree.  The  doubling  of  the 
errors  approximately  doubles  the  rmse's  of  all  the  estimated  variables. 

The  rmse  approximates  the  asymptotic  bias  of  the  estimator. 

Tables  4  and  5  present  the  simulation  results  when  the  range  to  target 
is  doubled  but  the  bearing  errors  are  0.2  and  0.4  degree.  With  the  excep¬ 
tion  of  the  estimated  range,  the  rmse's  of  the  estimates  are  approximately 
equal  to  those  obtained  when  the  errors  were  doubled.  Increasing  either 
the  range  to  target  or  the  errors  on  the  bearing  estimates  by  a  multipli¬ 
cative  constant,  k,  increases  the  bias  by  approximately  k.  This  empirical 
relationship  does  not  hold  for  the  estimated  range.  Although  the  rmse  of 
the  estimated  range  does  approximately  double  when  the  bearing  errors  are 
doubled,  it  increases  by  approximately  a  factor  of  4  when  the  range  to 
target  is  doubled. 

The  results  of  the  simulation  support  the  theory  presented  in  the 
paper.  These  results  are  a  function  of  the  particular  tracking  sequence 
chosen  as  well  as  the  target's  track.  The  circular  track  as  well  as  the 
90-degree  heading  were  chosen  for  the  simulations  because  of  their  simplicity. 

NOTE:  THESE  SYMBOLS  ARE  USED  IN  THE  TABLES  AND  ARE  IDENTIFIED  AGAIN 
FOR  CLARITY. 


rmse 

y.ji) 

(xB«  yB) 
N 

R 


root  mean  square  error  of . 

coordinates  of  the  target 
coordinates  of  the  tracker 
sample  size 

average  range  to  target 
rmse  of  e 

error  term  of  regression 
range  to  target 
heading  to  target 
velocity  of  target 

statistical  estimator  of  a  parameter 
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TABLE  3 


RESULTS  OBTAINED  USING  AD  HOC  ESTIMATOR  OF  RANGE  O  =  0.40 

£ 

Standard 


End  of  Track 

Actual 

Mean 

Mean  Bias 

Deviation 

RMSE 

xt(N/2) 

1.9799986 

1.9826841 

0.0026848 

0.2073205 

0.2073379 

yT(N/2) 

19.7999878 

19.9087372 

0.1087440 

0.1365973 

0.1745969 

Rt(N/2) 

19.8363953 

19.9413757 

0.1049771 

2.0171242 

2.0198536 

5t 

90.0000000 

89.8330841 

-0.1668751 

1.0594177 

1.0724792 

VT 

0.3600000 

0.3613726 

0.0013726 

0.0376254 

0.0376504 

Maximum  and  minimum  range  < 

during  tracking: 

20.181183 

and  19.485062 

N  =  45 

R  =  19.80 

rmseT(x 

T)  =  0.138 

rmseT(y1 

,)  =  1.334 

TABLE  4 

RESULTS  OBTAINED  USING  AD  HOC  ESTIMATOR  OF 


RANGE 

o  =  0.20  AND 
£ 

R  =  39.60 

Standard 

End  of  Track 

Actual 

Mean 

Mean  Bias 

Deviation 

RMSE 

xt(N/2) 

1.9799986 

1.9804792 

0.0004807 

0.2079040 

0.2079045 

yT(N/2) 

39.5999756 

39.7122803 

0.1123022 

0.1367894 

0.1769834 

Rt(N/2) 

39.6180878 

39.7253571 

0.1072704 

3.9831715 

3.9846144 

a 

90.0000000 

89.8725433 

-0.1274239 

1.1198845 

1.1271105 

VT 

0.3600000 

0.3608412 

0.0008411 

0.0373895 

0.0373989 

Maximum  and  minimum  range  during  tracking:  39.951111  and  39.281448 
N  =  45  rmse^,(xT)  =  0.138  rmseT(yT)  =  2.653 
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1 

1 


TABLE  5 

RESULTS  OBTAINED  USING  AD  HOC  ESTIMATOR  OF 


RANGE 

a£  =  0.40  AND 

R  =  39.60 

Standard 

End  of  Track 

Actual 

Mean 

Mean  Bias 

Deviation 

RMSE 

xt(N/2) 

1.9799986 

2.0484772 

0.0684777 

0.4705369 

0.4754934 

yT(N/2) 

39.5999756 

41.0848694 

1.4848890 

0.2057871 

1.4990807 

Rt(N/2) 

39.6180878 

41.1012878 

1.4831972 

9.1944504 

9.3133106 

<*T 

90.0000000 

89.9536896 

-0.0462814 

2.1883059 

2.1887951 

VT 

0.3600000 

0.3733910 

0.0133910 

0.0855601 

0.0866016 

Maximum  and  minimum  range  during  tracking:  39.951111  and  39.281448 


N  =  45  rmse^(x^)  =  0.276  rmse^Cy.j,)  =  5.306 
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